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The Transition from Anti-Parallel to Component Magnetic 
Reconnection 

M. Swisdak, 1 J. F. Drake, 2 M. A. Shay 2 and J. G. Mcllhargey, 3 

Abstract. We study the transition between anti-parallel and component collisionless 
magnetic reconnection with 2D particle-in-cell simulations. The primary finding is that 
a guide field wO.l times as strong as the asymptotic reconnecting field — roughly the 
field strength at which the electron Larmor radius is comparable to the width of the elec- 
tron current layer — is sufficient to magnetize the electrons in the vicinity of the x-line, 
thus causing significant changes to the structure of the electron dissipation region. This 
implies that great care should be exercised before concluding that magnetospheric re- 
connection is antiparallel. We also find that even for such weak guide fields strong inward- 
flowing electron beams form in the vicinity of the magnetic separatrices and Buneman- 
unstable distribution functions arise at the x-line itself. As in the calculations of Hesse 
et al. [2002] and Yin and Winske [2003], the non-gyrotropic elements of the electron pres- 
sure tensor play the dominant role in decoupling the electrons from the magnetic field 
at the x-line, regardless of the magnitude of the guide field and the associated strong 
variations in the pressure tensor's spatial structure. Despite these changes, and consis- 
tent with previous work, the reconnection rate does not vary appreciably with the strength 
of the guide field as it changes between and a value equal to the asymptotic reversed 
field. 



1. Introduction 

The fast dissipation of magnetic energy in collisionless 
plasmas is a common occurrence in nature, with examples 
ranging from tokamak sawtooth crashes to magnetospheric 
substorms to solar flares. The process common to these phe- 
nomena is thought to be magnetic reconnection, in which 
oppositely directed components of the magnetic field cross- 
link, forming an x-line configuration. The expansion of the 
newly connected field lines away from the x-line converts 
magnetic energy into kinetic energy and heat while pulling 
in new flux to sustain the process. 

Observations suggest that in many systems the ratio of 
the characteristic reconnection time to the Alfven crossing 
time is ~ 0.1. The simplest magnetohydrodynamic (MHD) 
description of reconnection [Sweet, 1958; Parker, 1957] is in- 
consistent with this value, being too slow by several orders of 
magnitude [Biskamp, 1986] . However, the numerical simula- 
tions comprising the GEM Reconnection Challenge [Birn et 
al., 2001] showed that the inclusion of the Hall effects, which 
are important at small spatial scales and are neglected in 
MHD, can produce fast reconnection. The magnetic topol- 
ogy in these simulations was understandably quite simple: 
equal and anti-parallel fields separated by a thin current 
layer. Yet even in the magnetotail, where this approxima- 
tion is often close to reality, a small field directed parallel to 
the current (a guide field) is often observed [Israelovich et 
al, 2001]. 

The effects of a guide field, B g , on magnetic reconnection 
have been examined before. Sharp differences have been 
seen in the large-scale flows around the x-line [Hoshino and 
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Nishida 1983; Tanaka 1995; Pritchett 2001] as well as the 
pressure and magnetic field signatures [Kleva et al, 1995; 
Rogers et al. 2003]. Three-dimensional particle simulations 
of similar systems without [Zeiler et al. 2002] and with 
[Drake et al., 2003] a guide field showed that the former 
was basically laminar in the direction parallel to the guide 
field while the latter developed strong turbulence. Pritch- 
ett and Coroniti [2004] noted that moderate guide fields 
(Bg/Bo < 1, where Bo is the reconnecting field) have only 
a slight effect on the reconnection rate, although Ricci et 
al. [2004] found somewhat slower rates for larger fields 
(B g /B = 3,5). Hesse et al. [1999,2002] and Yin and 
Winske [2003] showed that non-gyrotropic electron motions 
balance the reconnection electric field at the x-line in both 
the anti-parallel and guide field cases. 

In light of these results, determining the minimum guide 
field B g that changes the structure of the x-line becomes of 
interest. If it satisfies B g > Bo then the effects of a guide 
field can usually be ignored in the magnetosphere. On the 
other hand, if the transition occurs when B g « Bo guide- 
field reconnection is typical, and anti-parallel is a special 
case perhaps only relevant in simulations. We argue, based 
on both simulations and theoretical grounds, that the tran- 
sition occurs when the electron Larmor radius in the guide 
field at the x-line becomes smaller than the width of the 
electron current layer, i.e. for B g /Bo ~ 0.1. The implica- 
tion is that most magnetospheric reconnection is probably 
component reconnection. 

In section 2 of this paper we present our computational 
scheme and initial conditions. Section 3 presents results for 
the case B g — 0, section 4 for B g = 1.0 and section 5 for 



B 



0.2. We summarize our results and discuss their im- 



plications for understanding magnetospheric reconnection in 
section 6. 



2. Computational Details 

Our simulations are done with p3d, a massively parallel 
particle-in-cell code [Zeiler et al, 2002] . The electromag- 
netic fields are defined on gridpoints and advanced in time 
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with an explicit trapezoidal-leapfrog method using second- 
order spatial derivatives. The Lorentz equation of motion 
for each particle is evolved by a Boris algorithm where the 
velocity v is accelerated by E for half a timestep, rotated by 
B, and accelerated by E for the final half timestep. To en- 
sure that V • E = 4ivp a correction electric field is calculated 
by inverting Poisson's equation with a multigrid algorithm. 

The equations solved by the code are written in normal- 
ized units. Masses are normalized to the ion mass rm, the 
magnetic field to the asymptotic value of the reversed field, 
and the density to the approximate value at the center of the 
current sheet (see below) . Other normalizations derive from 
these: velocities to the Alfven speed va, lengths to the ion 
inertial length c/ui pi = di, times to the inverse ion cyclotron 
frequency fi"/, and temperatures to rmv A . 

Our coordinate system is chosen so that the inflow and 
outflow for an x-line are parallel to y and x, respectively. 
The guide magnetic field and reconnection electric field are 
parallel to z. For comparison, our x, y, and z unit vec- 
tors correspond to — x, z, and y in GSM coordinates. The 
simulations presented here are two-dimensional in the sense 
that out-of-plane derivatives are assumed to vanish, i.e., 
d/dz = 0. 

The initial equilibrium comprises two Harris current 
sheets [Harris 1962] superimposed on a ambient popula- 
tion of uniform density. The reconnection magnetic field is 
B x = tanh[(j/ — L y /4)/wo] — tanh[(?/ — 3L y /4)/wo] — 1, where 
Wo = 0.25 and L y — 6.4 are the half-width of the initial cur- 
rent sheets and the box size in the y direction respectively. 
This configuration has two current sheets and allows us to 
use fully periodic boundary conditions. The electron and 
ion temperatures, T e = 0.05 and Tj = 0.5, are initially uni- 
form as is the guide field B g . Except for the background 
(lobe) population, which can have arbitrary density (here 
n e — 0.2), pressure balance uniquely determines the ini- 
tial density profile. In this equilibrium the density at the 
center of each sheet is w 1.1 at t = 0. At t = we per- 
turb the magnetic field (B x /Bo ~ 0.1) to seed x-lines at 
(x,y) = (L x /4,3L y /4) and (3L B /4,L tf /4). 

To conserve computational resources, yet assure a suf- 
ficient separation of spatial and temporal scales, we take 
the electron mass to be 0.01 and the speed of light to be 
20. The domain measures 6.4 on a side and the grid has 
1024 x 1024 points, which implies that there are wl6 grid- 
points per electron inertial length and 2 per electron Debye 
length. To check for convergence we doubled the box size 
for one run (for B g — 0.2) and saw no significant variation 
in our results. 

The particle timestep is 6 x 10~ 4 , or 0.12cj pe . Our simu- 
lations follow ~ 10 9 particles and conserve energy to better 
than 1 part in 1000. 

3. Overview, B g = 

Investigating the critical value of B g with multiple 3-D 
simulations would entail a prohibitive computational ex- 
pense, so we instead performed a series of 2-D simula- 
tions that varied only in the strength of the guide field. 
The restricted dimensionality means that many turbulent 
modes, including the Buneman instability seen by Drake et 
al. [2003], are not present. However, other investigators 
have found that the gross morphological features of recon- 
nection x-lines are roughly invariant in the direction parallel 
to the current density [Pritchett and Coroniti, 2004]. 

A snapshot of the out-of-plane current density near the 
x-line for a simulation with B g — is shown in Figure 1. 
The initial current sheet has completely reconnected; the 
plasma at the x-line at the time shown was in the low den- 
sity (ne = 0.2) lobe at the simulation's beginning. Since 



|B| = at the x-line for B g = inward- flowing electrons 
must at some point find themselves in a region where the 
magnetic field is too weak to dominate their motion. This 
happens at a distance from the x-line given roughly by the 
electron inertial length, c/ui pe = d e , which for this simu- 
lation is y m e /ne w 0.2 (in normalized units). Once the 
electrons demagnetize they stream towards the x-line paral- 
lel to the y axis, through the region of low magnetic field, 
until they turn due to the increasing field on the opposite 
side of the current layer and reverse direction. They then 
execute "figure-8" trajectories [Speiser 1965], oscillating in 
the y direction until escaping from the ends of the layer. At 
the turning points of their trajectories (where v y — ► 0), the 
local electron density increases, forming a bifurcated cur- 
rent sheet. Zeiler et al. [2002] have previously reported 
this bifurcation, although it is particularly noticeable in our 
simulations because of the high spatial resolution (16 grid- 
points per d e ) and large ion to electron temperature ratio 
(Ti/T e = 10). The bifurcation would be obscured in simula- 
tions where these parameters had smaller values. Note that 
this bifurcation is at a much smaller scale than the ion-scale 
splits reported in Cluster observations [Runov et al. 2003]. 

The bifurcation is also evident in Figure 2, which shows 
the diagonal components of the electron temperature. In 
analogy with the definition of the fluid pressure tensor we 
define the electron temperature tensor in a grid cell as a sum 
over the N local particles: 

N 

T al 3 = <V a >){vp,i- <V/}>), (1) 

i— 1 

where <...> denotes an average, e.g., <x>= J2f=i x i- 
Like the pressure tensor, which is related to the temperature 
tensor by P a p — nT a p where n is the density, the temper- 
ature is symmetric, T a @ — Tp a . For an isotropic plasma 
the off-diagonal elements of the temperature tensor vanish 
while the diagonal elements are equal to each other and to 
the scalar temperature, T aa = T e . This is the case at t = 
in our simulations. 

The decrease in Tyy and T zz during inflow are consistent 
with the adiabatic invariance of the magnetic moment jj,: 
B « B x decreases, while fj, oc v±/B remains constant. T xx , 
approximately the parallel temperature, simultaneously in- 
creases due to energy conservation. Any energy change due 
to the interaction of the reconnection electric field E z and 
the curvature and grad-B drifts is small everywhere except 
near the x-line. 

Once inside the layer the electrons demagnetize and, as 
previously seen by Zeiler et al. [2002], the electron distri- 
bution in v y space separates into two counter-propagating 
beams due to the cross-current layer bounce motion (see 
Figure 2d). As a consequence T yy sharply increases, as has 
previously been noted by Horiuchi and Sato [1997]. 

Despite the beams we see no evidence of a two-stream in- 
stability, probably because of the small current layer width. 
Unstable wavenumbers for the electron-two-stream instabil- 
ity satisfy k y vo < 2u> pe where vq is the separation of the 
beam velocities [Krall and Trivelpiece 1986]. The maximum 
growth rate, 7 = uj pe /\/8, occurs for k y vo = \/3/2 uj pe , and 
as vo — > the growth rate vanishes. In the simulation the 
beam separation is largest at the x-line, vo « 8, and drops 
to at the edges of the layer w 0.22(= d e ) upstream. With 
a local Lu P e « 90 the instability criterion implies that only 
wavelengths A > 0.28 are unstable at the x-line, and thus 
that the two-stream instability is not excited in the narrow 
current layer. 

In a 2-D collisionless plasma the reconnection electric field 
at the x-line is ultimately balanced by the divergence of the 
electron pressure tensor [Vasyliunas 1975]. In our units the 
collisionless electron fluid momentum equation is 

E = -v e xB - — V-P e - m e (v e -V)v e - (2) 

n e at 
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which is exact insofar as the pressure tensor P e incorpo- 
rates all kinetic effects not included in the other terms. The 
reconncction electric field is thus 



-(V X By — VyB X ) 



dP x: 
dx 



+ 



dP v . 



dy 



(3) 



dv z 



dv z dv z 



where we have dropped the electron subscript and used the 
fact that d/dz = 0. In a steady state only the pressure terms 
can balance E z at the x-line, as can be seen in Figure 3. Far 
from the current layer the EMHD relation E = — v e x B 
holds, while nearer the x-line both the off-diagonal elements 
of the pressure tensor and the inertial terms are important. 
At the x-line the pressure tensor terms dominate. Both 
dPy Z /dy and dP xz /dx contribute, although the former is 
larger in our simulation by a factor of w 2. The term pro- 
portional to d/dt is not shown separately but, as expected 
during quasi-steady reconnection, is negligible. 



Conservation of magnetic moment again explains the de- 
crease in T± along the inflow direction. Since B g is constant 
the decrease in B is smaller than was the case in section 3, 
and the relative decrease of T± in Figure 6c is smaller than in 
Figure 2d. Because the electrons remain magnetized at the 
x-line, T± remains small, in sharp contrast to the results 
shown in Figure 2. The increase in Tj| inside the current 
sheet is due to the intermixing of colder inflowing electrons 
and electrons accelerated by the parallel electric field along 
the separatrices. The unimodal distribution function of Fig- 
ure 6d confirms that the electrons do not execute Speiser-like 
orbits. 

The terms balancing the reconnection field in equation 
(3) are shown for this case in Figure 7. At the x-line the 
off-diagonal elements of the pressure tensor again make the 
primary contribution although, as a comparison of Figures 
7 and 3 makes clear, the scale length over which they are 
important is much smaller than when B g = 0. This is con- 
sistent with the results of Hesse et al. [2002]. Unlike the 
anti-parallel case the dP xz /dx term makes the dominant 
contribution while the dP yz /dy term is negligible. 



4. Guide Field, B g = 1.0 

A large guide field changes the structure of the x-line by 
both lowering the total plasma f3 and magnetizing the elec- 
trons throughout the current sheet. When B g is small the 
dominant wave mode at small lengthscales, and hence the 
governor of the particle dynamics, is the whistler. This is 
seen when electrons in the outflow region are accelerated 
by E z and drag the magnetic field out of the reconnection 
plane [Mandt et al., 1994; Shay et al., 1998], causing the 
well-known quadrupolar symmetry in B z along the separa- 
trices [Sonnerup, 1979; Terasawa, 1983] . As B g increases 
the importance of the kinetic Alfven mode grows [Rogers 
et al, 2001]. For that mode the coupling occurs when E\\ 
accelerates electrons along newly reconnected field lines, in- 
creasing the electron density on one side of the current layer, 
decreasing it on the other, and forming a quadrupolar pat- 
tern [Kleva et al, 1995]. The perturbations in B z acquire 
a component determined by pressure balance, leading to 
a symmetric component that can, for very large B g , over- 
whelm the quadrupolar pattern [Rogers et al, 2003]. Figure 
4 shows the electron density and B z from the simulation 
discussed in the previous section (B g = 0) and one that is 
otherwise identical except that B g — 1.0. 

The parallel velocity of the electrons, mostly directed out 
of the reconnection plane, develops a quadrupolar symmetry 
opposite to that of the density (high density paired with low 
velocity and vice versa). The density asymmetry is a larger 
effect, however, and the result is an out-of-plane current den- 
sity that is canted with respect to the initial current sheet, 
as can be seen in Fig. 5. Because inflowing electrons remain 
magnetized in the guide field they do not have figure-8 tra- 
jectories at the x-line and the bifurcations in the electron 
density and current density disappear. 

In the B g — simulation discussed in Section 3 the mag- 
netic field on the inflow axis (x — 4.8) was dominantly par- 
allel to x (except at the x-linc where |B| = 0). For B g = 1.0, 
in contrast, the field rotates « 45° between the lobe plasma 
and the x-line. The rotation complicates the interpretation 
of the temperature tensor of equation (1), so to simplify we 
transform to a coordinate system where the axes are parallel 
and perpendicular to the local magnetic field. In this frame 



T= [T N bb + Tx(I-bb)] + T ng 



(4) 



where b is a unit vector in the direction of the magnetic field, 
I is the unit tensor, and T ng contains the non-gyrotropic 
terms. Images of the parallel and perpendicular temper- 
atures along with cuts through the x-line are shown for 
B g = 1.0 in Figure 6a-c. 



5. Transition 



tion of Section 3 and that of Section 4 occur? Far from the 
x-line, where both species are completely magnetized, the 
guide field cannot play an important dynamical role unless 
Bg > Bo- As a particle approaches the current layer, how- 
ever, the reconnecting component decreases and the influ- 
ence of the guide field rises. Qualitatively, B g is important 
when the associated electron Larmor radius is equal to the 
spatial scale associated with the x-line. 

Consider a system with B g = and examine the z compo- 
nent of the electron equation of motion under quasi-steady 
conditions, 



dv ez 

m e v ey — — = -eE z 
dy 



-(vexBy - v ey B x ) - — ®^ yz , (5) 
c n dy 



where we have assumed that derivatives with respect to x 
can be neglected when compared to those with respect to 
y. In a 2D steady-state system Faraday's Law implies that 
E z is relatively uniform (see Figures 3 and 7). Far from 
the x-line electrons are frozen to the magnetic field and the 
E z and v x B terms are roughly equal. Within the current 
layer the convective part of the inertial term and the pres- 
sure tensor become important, with the transition occurring 
at some lengthscale A y where the terms balance. If, for sim- 
plicity, we restrict our attention to the vertical axis through 
the x-line, symmetry implies that B y is zero and 



(6) 



where (] IiUp = eB X}Up /m e c is the cyclotron frequency based 
on the reconnecting field at y — ±A y . In general B XiUp will 
be less than the asymptotic reconnecting field. Within this 
inner scale the electrons carry most of the current, so we 
also have 



47T 4tt 
— J z — — nv e 

c c 



dB 
dy 



V x B w — - (7) 



where we have ignored both the displacement current and 
the contribution from the current due to 8B y /dx. Con- 
verting derivatives with respect to y to division by A y and 
combining equations (6) and (7) we find that 



Ay = 



(8) 
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The electron velocity producing the current is v ez ~ d e Q x , up . 

Now consider the addition of a small ambient guide field. 
The magnetic field no longer vanishes at the x-line and the 
electron Larmor radius there is 

Po = 0) 

where we have taken v ex w based on symmetry consider- 
ations, f2 e , g = eBg/m e c is the electron cyclotron frequency 
based on B g , and v ey is the electron inflow velocity into the 
unmagnetized region around the x-line. The guide field will 
be important when this Larmor radius is smaller than the 
width of the current layer, p g < A y . 

The major contributors to the electron inflow velocity v ey 
are the thermal speed and the E x B drift. A Sweet-Parker- 
like scaling suggests that the ratio of the E x B inflow speed 
to the outflow speed is equal to the normalized reconnec- 
tion electric field E z . For a wide range of conditions it has 
been shown [Shay et al, 1999; Shay et al, 2001] that the 
outflow is roughly equal to the electron Alfven speed vae 
and E z w 0.1. In the low-temperature limit (vth « vexb) 
this argument implies that v ey w O.lvAe and the bound for 
a dynamically important B g is given by 

V ey 0.1«Ae \ , „ n , 

P9 = tt- = -o < A y = d e (10) 

or 

B g >E z ~ 0.1. (11) 

In our simulations E z , and hence the transitional value of 
B g , is w 0.2. 

Equation (11) is not valid when the electron thermal 
speed v t he dominates the contribution from the E x B drift. 
In the high-temperature limit one must substitute v t h rather 
than vexb for v ey in equation (9), and the relevant crite- 
rion becomes f} g < 1 where f3 is evaluated with the lobe 
density and temperature. In the simulations presented here 
Vthe ~ vexb and so the critical value of the guide field re- 
mains B g « 0.2. 

We note that the estimate for the electron inflow velocity 
v ye w O.Icab is smaller than the counter-streaming velocity 
of the electrons shown in Fig. 2(d). This is because a local 
electrostatic field E y develops inside the electron current 
layer that accelerates the electrons towards the magnetic 
null. However, this field decelerates the electrons once they 
cross the null so this electric field does not change the ef- 
fective electron Larmor radius of the electrons in the guide 
field. 

We have explored the transition from anti-parallel to fi- 
nite guide field reconnection through a series of simulations 
with Bg = 0, 0.1, 0.2, 0.4, and 0.75 and found, in agreement 
with our above arguments, that simulations with B g = 0.2 
most clearly display characteristics intermediate between 
B g = and B g = 1. 

Figure 8 shows the out-of-plane current density for a run 
identical to those previously discussed except that B g — 0.2. 
There is no bifurcation in the current density and the canting 
of the current layer, while present, is not as strong as the case 
with B g = 1.0 (Figure 5). Figure 9 shows the electron den- 
sity and B z for this simulation. The electron density is not 
bifurcated and bears some resemblance to the quadrupolar 
cavities so prominent in Figure 4c, while B z , although still 
basically quadrupolar, no longer has the strong symmetry 
obvious of Figure 4b. Evidence for a transition can also be 
seen in the parallel and perpendicular temperatures shown 
in Figure 10. The results are clearly intermediate between 
Figures 2 and 6. The cut in Figure 10c shows that the in- 
crease in the parallel temperature as electrons approach the 
x-line is similar in both magnitude and profile to the B g = 



case. Within the current layer, however, Tji increases to a 
sharp peak similar to that for B g = 1. The perpendicular 
temperature decreases towards the x-line and then, within 
the layer, rises to a peak. This peak is midway in magnitude 
between the B g = and B g = 1 cases. 

In order to examine the variation with guide field of the 
off-diagonal pressure tensor terms of Equation 3 it is neces- 
sary to separate the gyrotropic and non-gyrotropic contri- 
butions. The gyrotropic part is strongly influenced by the 
presence of a guide field and, in any case, does not contribute 
to balancing E z at the x-line. Figure 11 shows the change in 
the non-gyrotropic portion of P yz as the guide field varies. 
The cuts in Figure lid demonstrate that as B g increases the 
role of dP yz /dy in balancing the reconnection electric field 
at the x-line decreases dramatically. 

The presence of a small guide field also has signatures far 
from the x-line. Figure 12 shows the v x -v y distribution func- 
tions for three simulations with different guide fields taken 
at the same point, just upstream of the upper-left separatrix 
and, for the runs with a finite guide field, inside the density 
cavity. Because of the low plasma density within the cav- 
ity the parallel electric field remains finite over an extended 
region along the separatrix [Pritchett and Coroniti, 2004]) 
in the B g = 1.0 system. This electric field locally acceler- 
ates the electrons, producing a strong beam flowing towards 
the x-line. However even for B g — 0.2 the beam is already 
clearly present. (The features in the upper left quadrant of 
each panel are electrons that have already been accelerated 
at the x-line). There is a net current but no distinct beam 
for B g = 0. The origin, detailed structure, and effect of this 
extended region of E\ \ is discussed in a future publication. 

6. Discussion 

This study suggests that only a minimal guide field, 
Bg « O.lBo is required to alter the dynamics of electrons 
both in the vicinity of the x-line and at remote locations 
along the separatrices. The implication is that in most real 
systems, including the magnetotail, the guide field might 
not be negligible. In any case this study suggests that one 
can not simply ignore the guide field if B g <JC Bo. 

A counter-argument could be made that reconnection 
with B g 7^ is significantly slower compared to the case 
with Bg — and therefore the magnetosphere will self-select 
locations where the guide field is nearly zero (less than 0.1 
of the anti-parallel field). We find that this is not the case. 
The guide field alters the dynamics both locally and at large 
scales significantly before the rate of reconnection is signifi- 
cantly affected. Specifically, Figure 13 shows that magnetic 
flux reconnects only slightly (« 10%) slower for B g = 1.0 
than for B g — 0, a result consistent with other simulations 
[Rogers et al, 2003; Pritchett and Coroniti, 2004]. It should 
be noted, however, the the onset of reconnection in real sys- 
tems could be biased either for or against guide fields. Our 
simulations do not address this question since they start 
with a finite perturbation that effectively places the system 
in the nonlinear regime at t — 0. 

Another factor potentially affecting guide field reconnec- 
tion is the effect of an ambient pressure gradient and the 
associated diamagnetic drifts. At the magnetopause, where 
density gradients perpendicular to the current layer produce 
diamagnetic drifts, the reconnection rate can be strongly re- 
duced [Swisdak et al., 2003]. However, it was shown that dia- 
magnetic suppression occurs for small guide fields with the 
transition occurring when (3 X > B g /Bo (for density length 
scales of order an ion inertial length) . Combined with the re- 
sults of this work the conclusion is again that magnetopause 
reconnection always includes a dynamically important guide 
field. 

Guide fields also play an important role in the develop- 
ment of turbulence in three-dimensional reconnection simu- 
lations. Simulations by Drake et al. [2003] with B g — 5.0 
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showed that turbulence can self-consistently develop at a 
reconnection x-line. The acceleration of electrons by the re- 
connection electric field led to a separation of the ion and 
electron drift speeds, which then triggered the Buneman in- 
stability. At late time the nonlinear evolution led to the for- 
mation of electron holes, localized bipolar regions of electric 
field. These structures produced an effective drag between 
the ions and electrons that was large enough to compete 
with the off-diagonal pressure tensor in balancing the recon- 
nection electric field. Earlier 3D simulations with B g = by 
Zeiler at al. [2002] produced no significant turbulence at the 
x-line once reconnection was established. The strong elec- 
tron heating for B g — (as shown in Figure 2) suppressed 
all streaming instabilities near the x-line since for T; 3> T e 
such instabilities require a beam velocity vt greater than 
the electron thermal speed v te - Yet since B g = 5.0 is sig- 
nificantly larger than typical magnetospheric values it was 
unclear from these studies whether turbulence and enhanced 
ion-electron drag were common features of reconnection in 
the magnetosphere. 

Our 2D simulations cannot produce the Buneman insta- 
bility at the x-line seen in these earlier simulations. How- 
ever, we can examine the distribution functions produced 
by our simulations and determine whether they would be 
unstable in a full 3D system. In the low temperature limit 
with k parallel to both B and the relative drift velocity Vb 
a plasma is Buneman unstable to wavenumbers satisfying 
the relation kvb < oJ pe [Krall and Trivelpiece, 1986]. For fi- 
nite temperature plasmas with Maxwellian distributions the 
condition for instability is more complicated and, in fact, for 
small drifts in warm plasmas no instability exists. The in- 
stability threshold for arbitrary distributions can be found 
numerically, but a rough rule of thumb is that a plasma is 
Buneman unstable if the electron and ion velocity distribu- 
tion functions do not substantially overlap (vt, > v te ). 

Figure 14 shows the distribution of v z at the x-line for 
runs with three values of the guide field, B g = 0,0.2,1. 
The dotted line shows the ion distribution and the double 
peak is a remnant of the two initial populations, the non- 
drifting background and the drifting population comprising 
the initial current sheet. The three electron distribution 
functions demonstrate that as the guide field increases the 
x-line electrons acquire a larger drift with respect to the 
ions. The electron energy gain is limited by the amount of 
time they spend within the current layer before advecting 
into the outflow region. Larger simulations suggest that our 
small system size may prevent the electrons from reaching 
their maximum speeds. A finite B g acts like a guide wire, 
restraining this flow and leaving more time for a particle 
to be accelerated by the reconnection electric field. If our 
simulation were three-dimensional the B g — 1 case would 
almost certainly be Buneman unstable, while the B g = 
case would not. Again, the B g — 0.2 case is transitional. 

The small size of the transition guide field has important 
implications for magnetospheric reconnection. The x-line 
current sheet should usually be canted with respect to the 
ambient current layer. Density and flow velocities measured 
on the separatrices should have a quadrupolar symmetry. 
Distribution functions taken just upstream of the separatri- 
ces should exhibit an inward- flowing beam. Turbulence, and 
electron holes in particular, should be common. Magneto- 
spheric reconnection with negligible guide fields should be 
rare. 
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Figure 1. Reconnection in a system with B g — at 
t = 4.5. (a) The out-of-plane current density overlaid 
with magnetic field lines in a region surrounding the x- 
line. The blackest regions are not color-coded correctly, 
having been over-exposed to show details near the x-line 
(b) A vertical cut through the x-line at x = 4.8. 
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Figure 2. Data for B g — 0. Panels (a) and (b) show 
the electron temperatures T xx and T yy (see the text for 
definitions) near the x-line of Figure 1. Panel (c) shows 
cuts of the temperatures at x — 4.8. The solid line is T yy , 
the dashed T xx , and the dotted T zz . Panel (d) shows the 
distribution of v ey in a region measuring 1.0 d e long and 
0.5 d e high and centered on the x-line. 
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Figure 3. For B g = a cut through the x-line at x = 
4.8 showing the various terms balancing E z in equation 
(3). E z is shown by the stars and the dashed, solid, and 
dotted lines denote the v x B, divergence of the pressure 
tensor, and inertial terms, respectively. To reduce the 
noise the plotted quantities were averaged over 20 grid- 
points in x and 4 in y. 






Figure 4. The electron density and out-of-plane mag- 
netic field B z in simulations with B g — and 1.0. (a) 
ne, B g = 0; (b) B z , B a = 0; (c) n e , B g = 1.0; (d) B z , 
B g = 1.0. In all panels t = 4.5. In panels (a) and (c) 
black areas have been over-exposed to show detail at the 
x-line. 
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Figure 5. Reconnection in a system with B g — 1.0 at 
t = 4.5. (a) The out-of-planc current density overlaid 
with magnetic field lines in a region surrounding the x- 
line. Black areas have been over-exposed to show detail 
at the x-line. (b) A vertical cut through the x-line at 
x = 4.8. 
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Figure 8. Reconnection in a system with B g — 0.2 at 
t — 4.5. (a) The out-of-planc current density overlaid 
with magnetic field lines in a region surrounding the x- 
line. Black areas have been over-exposed to show detail 
at the x-line. (b) A vertical cut through the x-line at 
x = 4.8. 




Figure 6. Data for B g = 1.0. Panels (a) and (b) show 
Tj| and T±, respectively. Panel (c) shows vertical cuts at 
x = 4.8. The solid line is Tji and the dashed T±. Panel 
(d) shows the distribution of v ey in a region measuring 
1.0 d e long and 0.5 d e high and centered on the x-line. 
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Figure 9. The electron density (a) and out-of-plane 
magnetic field B z (b) in a simulation with B g = 0.2. 
In panel (a) black areas have been over-exposed to show 
detail at the x-line. 



Figure 7. The terms balancing the reconnection elec- 



dashed, solid, and dotted lines denote the v x B, diver- 
gence of the pressure tensor, and inertial terms, respec- 
tively. Note the difference in horizontal scale between 
this figure and Figure 3. To reduce noise the plotted 
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Figure 10. Data from B g = 0.2. Panels (a) and (b) 
show Tj and T±, respectively. Panel (c) shows vertical 
cuts at x = 4.8. The solid line is Tji and the dashed T±. 
Panel (d) shows the distribution of v ey in a region mea- 
suring 1.0 d e long and 0.5 d e high and centered on the 
x-line. 
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Figure 11. Panels (a)-(c) show the non-gyrotropic com- 
ponent of P yz for B g = 0,0.2, and 1.0. Panel (d) shows 
a cut at x = 4.8 in each figure. The solid, dashed, and 
dotted lines are for, respectively, B g — 0,0.2, and 1.0. 





Figure 12. Panels (a)-(c) show the v x -v y distribution 
functions taken just upstream of the separatrix for runs 
with B g = 0,0.2,1.0, respectively. The slanted line in 
each panel shows the direction of the local magnetic field. 
The colors are measured in units of phase space den- 
sity. The counter-propagating features at the upper left 
of each figure are due to electrons accelerated at the x-line 



1.0 

o.s 

[Z 

ti 0.6 
m 

I 

I 0-4 
0.2 
0.0 



ii 



Figure 13. The reconnected flux versus time for the 
three runs discussed in this work. The solid, dashed, and 
dotted lines correspond to B g = 0, 0.2, and 1.0, respec- 
tively. The reconnection rate is the slope of the curves. 




Figure 14. The v z distribution functions for three runs. 
The solid line i 
B g = 0.2 and 
shows the ions. 



B g = 0.2 and the dashed B g = 1.0. The dotted line 



